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ABSTRACT 



Eisenstein Loeb and Turner (1996) (ELT) have recently proposed a method for 
estimating the dynamical masses of large-scale filaments, whereby the filament is 
modeled by an infinite, axisymmetric, isothermal, self-gravitating, radially virialized 
cylinder, for which ELT derive a global relation between the (constant) velocity 
dispersion and the total line density. We first show that the model assumptions of ELT 
can be relaxed materially: an exact relation between the root-mean-square velocity and 
the line density can be derived for any infinite cylinder (not necessarily axisymmetric), 
with an arbitrary constituent distribution function (so isothermality need not be 
assumed). We then consider the same problem in the context of the modified dynamics 
(MOND). After a brief comparison between the scaling properties in the two theories, 
we study two idealized MOND model filaments, one with assumptions similar to those 
of ELT, which we can only solve numerically, and another, which we solve in closed 
form. A preliminary application to the same segment of the Perseus-Pisces filament 
treated by ELT, gives MOND M/L estimates of order 10(M/L) Q , compared with the 
Newtonian value M/L ~ 450(H /100kms- 1 Mpc~ 1 )(M/L) Q , which ELT find. In spite 
of the large uncertainties still besetting the analysis, this instance of MOND application 
is of particular interest because: 1. Objects of this geometry have not been dealt with 
before. 2. It pertains to large-scale structure. 3. The typical accelerations involved are 
the lowest so far encountered in a semi-virialized system-only a few percents of the 
critical MOND acceleration-leading to a large predicted mass discrepancy. 



Subject headings: Dark matter, modified dynamics, large-scale structure, superclusters, 
the Perseus-Pisces supercluster 



1. Introduction 



In a recent paper Eisenstein Loeb and Turner (1996)| (ELT) describe a method for estimating 
the dynamical masses of large-scale filaments. Elementary considerations lead one to expect that 
when such a filament is virialized in the radial direction its total line density (mass-per-unit-length) , 
fio, and its typical velocity dispersion, a, are related by G/j,q ~ a 2 . For concreteness, ELT consider 
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a model of these filaments consisting of an infinitely long, axially symmetric, self-gravitating 
cylinder whose constituent velocity dispersions are constant (we call it isothermal, even if the 
velocity distribution need not be "thermal"). For such cylinders, ELT derive the exact relation 

G^o = 2ai, (1) 

where a± is the velocity dispersion along a line of sight perpendicular to the cylinder axis, 
averaged over the whole cylinder. ELT study carefully the applicability of this relation by testing 
its performance on N-body-simulated filaments. Finally, they apply this relation to a segment of 
the Perseus-Pisce filament and estimate for it a B-band M/L ratio of about 45O(M/L) , implying 
large quantities of dark matter. 

The modified dynamics (MOND) was proposed as an alternative theory to Newtonian 



dynamics with a view to eliminating the need for dark matter [ Milgrom (1983) |. It posits that in 
the limit where the typical accelerations in a system governed by gravity are much smaller than 
some value (a ~ 10 -8 cm s~ 2 )-the deep-MOND limit-the mass, M, radius, r, and acceleration, 
a, are related by a 2 /a a ~ MG/r 2 , and not by the standard a ~ MG/r 2 . So, the smaller the 
mean acceleration in a system, the larger the expected mass discrepancy deduced by a Newtonian 
analysis. 



MOND has been tested extensively on galaxies [Begeman, Broeils, and Sanders (1991 



Sanders (1996) |, and it is natural to check how well it does in the context of large-scale structure. 



As explained in Milgrom (1989)| , MOND is not developed enough to afford application to most 



large-scale structures, with over-density not much exceeding unity, or where the Hubble flow is 
important. Because large-scale filaments seem to be at least radially virialized, they offer a unique 
opportunity to apply MOND to large-scale structure. Inasmuch as a filament is not completely 
virialized, and extends to radii where the Hubble flow is important, our analysis can be valid only 
at smaller radii. 

Large-scale structures constitute an extreme from the point of view of MOND because they 
involve accelerations that are the lowest in the range so far observed, and they carry the usual 
importance and interest of extremes: Values as low as 0.1a o are measured at the outskirts of dwarf 
spirals [e.g. Banders (1996) |, while, as we shall see below, in the Perseus-Pisces filament a is only 



a few percents of a Q . 

Section 2 brings some generalities concerning cylindrical and spherical structures in Newtonian 
dynamics and in MOND, and extends the ELT result to more general cylinders. In section 3 we 
describe two models for self-gravitating infinite cylinders in MOND. In section 4 we consider the 
Perseus-Pisces filament. 



2. Virial relations for cylindrical and spherical systems in Newtonian dynamics and 

in MOND 
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Virialized spherical systems of size r and velocity dispersion a, have a mean acceleration, 
~ cr 2 /r, that roughly equals MG/r 2 in Newtonian dynamics, and (MGa /r 2 ) l l 2 in MOND 
(very-low-acceleration limit). Thus, in the former the virial relation M ~ a 2 r/G contains the 
system's radius, while in the latter the analogous relation M ~ a 4 /Ga a does not. For cylinders, 
the tables are turned, and since the Newtonian mean acceleration is ~ fiG/r, while in MOND it is 
~ (nGdo/r) 1 / 2 , we have in Newtonian dynamics [iG ~ a 2 , while in MOND we have [iG ~ o" 4 /ra a . 

We see then that the Newtonian cylindrical case is akin to the spherical, deep-MOND case, 
the crucial common property being the logarithmic behavior of the potential at infinity. In fact, 
the isothermal axisymmetric cylinder discussed by ELT, and the deep-MOND, isothermal spheres 



described in Milgrom (1984) are special cases, with dimensions two, and three, respectively, of 
D-dimensional, self-gravitating, isothermal spheres, held together by a potential, cp, that satisfies 
a Poisson-like equation 

V-[\V<p\ a V<p]=Ap, (2) 

where a = D — 2 is chosen such that the potential of an isolated mass is logarithmic at infinity 
(for D = 2 we get the usual Poisson equation, for D = 3 we get deep MOND). The general density 



law of such spheres is ( Milgrom (1996 



(r/rn)-^- 1 ^ 

* r — a ,D (3) 

l + (r/r )D-i P 

{(5 is the anisotropy parameter-see the next section). These spheres satisfy a mass- velocity- 
dispersion relation of the form 

(vY' 1 = MAa' 1 (^Y'\ (4) 

where (v 2 ) is the mass-weighted root-mean-square velocity of the system, a D is the D-dimensional 
solid angle (2ir for D = 2, etc.). In fact, Milgrom (1996)| proves that this relation is not limited 



to isothermal spheres, but holds for an arbitrary virialized system subject to the above potential 
equation. 

This relation for deep-MOND, isothermal spheres {D = 3, A = AnxGa Q , 03 = 47r): 

MGa = ^(v 2 ) 2 , (5) 



was first derived in Milgrom (1984); it was later generalized, first by perhard and Spergel (1992) to 



spherically symmetric systems with general velocity distributions, and then, further, by |Milgrom 



(1994) to arbitrary (not necessarily spherical) systems. 



Because the results of Milgrom (1996) are not published yet we now give the proof of 



relation(|4|), for the (relatively simple) special case D = 2, which constitutes a generalization of 



the ELT result to arbitrary cylindrical systems. Following Milgrom (1994) mutatis mutandis, we 
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consider a stationary, self-gravitating system, symmetric under z translations, that is composed 
of various particle species with masses and distribution functions f k (f,v), f, and v being the 
position, and velocity in the x — y plane. As usual, take the second time derivative of the quantity 

Q = g J ^ r d2y mk r2 > ( 6 ) 

k 

which must vanish due to the stationarity of the system (if the system is not stationary but still 
bound its long-time average vanishes): 

£ / i\ i\ ,n t f k (f, „-) *r m t Mr, 1, 0, (7) 

k k 

The first term in eq.(0) is the mass- weighted, 2-D, mean-square velocity, (v 2 ), multiplied by 
/iQ-the total line density of the cylinder. In the second term we put a = —V(p, where <p is the 
gravitational potential, to obtain 

A*o(u 2 ) = Y d2r d2y m k fh(r,v) r- V<p. (8) 
k J 

The v integration, and the sum over species can now be performed to yield the standard result 



H (v A ) = J d 2 r p(r)r- V<p. (9) 
Now, using the Poisson equation we can write eq.(||) as 

Vo(v 2 ) = 4^ / ■ V<^V, (10) 

which can be written as 

Mv 2 ) = J V • [r ■ - i(V^) 2 r1d 2 r. (11) 

Applying Gauss's theorem to write the integral as a surface integral at infinity, and remembering 
that in the relevant geometry V</? goes asymptotically as — ► —2Gfiof/r 2 , the integration gives 
G/Uq, and we finally get the general virial relation 

{v 2 ) = G» . (12) 

When the system is not stationary, but still bound, (v 2 ) in eq.(jl2|) is replaced by its long-time 
average, {v 2 }. 

Note that, in fact, in deriving eq.(|i~2|) we do not make use of the perfect cylindricity of the 
system. We only use the Poisson equation, and the cylindrical boundary behavior of the potential 
at infinity. All still hold if the object looks like a uniform, straight line at infinity, i.e., it can 
have nonuniformities of different kinds: wiggles, modulations of the line density, etc., as long as 
it is uniform on average. The problem is that in this case the system cannot be expected to be 
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stationary, so we cannot use the instantaneous, observed value of (v 2 ), and we cannot easily relate 
it to a line-of-sight velocity dispersion, as we do below. 

The 2-D root-mean-square velocity, (v 2 ), of an astronomical system cannot be measured 
from our vantage point, and we would like to express it in terms of the one-dimensional, 
system-integrated, root-mean-square of the velocity component along a line of sight perpendicular 
to the cylinder axis, a±. This can be done, for example, when a± is independent of the azimuthal 
viewing angle, in which case 

4 = \(v 2 ), (13) 

and we may then write 

2a\ = Gfi , (14) 



which is the relation derived by ELT for their model filament. Two instances in which eq,(|13D is 
valid are: a. The velocity ellipsoid is isotropic everywhere (e.g. in a gaseous system), but the 
system need be neither isothermal, nor axisymmetric. b. The system is axisymmetric, but then 
the velocity distribution need be neither isotropic, nor isothermal. The ELT model is a special 
case of instance b. 

Another point to be made about the ELT model is that it can be extended, at no further 
cost, to include a very dense core at the center of the filament. If we approximate this core by a 
line singularity of finite line density /x(0), then the structure equation of ELT is modified only by 
replacing the anisotropy parameter (3 (see next section or ELT for denotation) by (3 = (3 + m,Q, 
where = 2p(0)G/a 2 . The behavior of the density near the origin is now p cx . convergence 
of the mass both near the origin, and at infinity is assured for (3 < 2, which sets an upper limit 
on the singular line density: /i(0) < (2 — f3)a 2 /2G. The relation between the "gas" line density 
(excluding the singularity) and the "gas" velocity dispersion is now 

2a 2 

ft (H = ^-2M0). (15) 

and the total line density is 

Hot = - M(0). (16) 
Our general relation(14) is seen to hold when we note that we must take in it the mass- weighted 



mean of <r^, including the contribution of the line singularity. This latter must itself have a 
a\(line) = fi(0)G/2. The line singularity can stand, for example, for an ELT model with a radius 
scale much smaller than that of the "envelope" gas. Such an extended model can help study 
possible departures from the ELT model. For example, when the "gas" mass is negligible the 
singularity line density determines a± through fi(0) = o\jG. 



3. Two idealized model filaments in MOND 
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Obtaining a quantitative inkling of the MOND relation between line density, radius, and 
velocity requires specific models. We now consider two such filament models. Not enough is known 
about the velocity distribution, and radial-density structure of real filaments to assess the degree 
of their adequacy. Even if in their very outer parts actual filaments depart greatly from these 
models, these may be applicable in the inner parts. 

The first model makes the same assumptions as ELT's. It consists of an infinitely long, 
axisymmetric, self-gravitating system such that the azimuthal, and radial velocity dispersions-^, 
and a r respectively-are position independent (accelerations along the axis are assumed negligible). 
We assume that all accelerations are much smaller than a Q so that the deep MOND limit is taken. 
The MOND Jeans equation can then be written as 



din p 3 
ar r 



din p 3 
dr r 



/i(r)G 



(17) 



Here, p(r) is the mass density at cylindrical radius r, p,(r) is the line density within radius r: 



p(r) = 2-7T / p(r')r' dr' , 



and 



/3 = l-o$/o* 



(18) 



(19) 



is the anisotropy parameter. The system-integrated rms velocity component along a line of sight 
perpendicular to the cylinder axis, a±, is given by 



4 



(20) 



Equation(|j"7|) is derived along the lines detailed in Milgrom (1984), and is based e.g. on the 
formulation of MOND in |Bekenstein and Milgrom (1984) . 



We shall see that the solution of eqs.(|T7|)(|T8|) leads to a finite value of the asymptotic p. 
Define then po = p(oo), and use it to define length, and density scales by 



a,. 



ro 



a o p G 



Po 



Po 
2irr. 



2 • 



respectively. The deep-MOND Jeans equation can then be written as 

c ,, -«r 1 c / +»r 1/2 c 1/2 c , = o, 

where s = 1 — (3, and we use the dimensionless variables 

n = r/r , X = p/p , 

and 

C(n) = p(nr )/p = / A(r/)?/ drj . 



(21) 

(22) 

(23) 
(24) 
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Equation(f22j) is to be solved with the boundary conditions £(0) = 0, and C(°°) = 1- 



Near the origin, the third term in eq.(22) can be neglected (it is the only term by which the 
Newtonian equation differs from the MOND equation), and we can solve the remaining equation 
to obtain there 

C(r « 0) « brf, ((r « 0) « -4rrrf +1 . (25) 

s + 1 

Near the origin thus p(r) oc r~@, as in the Newtonian case. Values [3 < give nonmonotonic 
density laws, and we ignore these. (The distribution function for such cylinders can be written as 
a function of the particle angular momentum, J = rat, and the energy, E = v 2 /2 + tp(r), for unit 
mass particles, f(v,r) oc J~^exp(— E/2a 2 ), and is increasing with J for (3 < 0.) 

At large rj the second term is negligible, and we may put £ ~ 1 in the third term. Equation (p2|) 
than gives for the asymptotic form of the density X(rj — ► oo) oc r/ _1 exp(— 2r/ 1/2 ), or 

p oc r- 1 exp[-2(r/r ) 1/2 ]. (26) 



From eq.(p5|) we see that £ is increasing near the origin. By eq.(|22|) it cannot decrease at 
larger radii because if it has a maximum and hence Q" vanish there and £ remains constant 
from that point hence; it can also not increase indefinitely as can be seen from eq.(^), and must 
then go to a constant at infinity, has we have assumed all along. 



The solution for £ can be found, numerically, by integrating eq.(|22|) out from the origin, 
starting with expressionpB]), shooting to find the value of the parameter b for which C(°°) = 1- 
We can avoid this trial- and-error if we rescale the variables rj, and £: 

rj = jfj, C = wC, (27) 

such that a. our structure equation (22) remains of the same form in fj, and and b. near the 



origin dt^/df) = t) s . This can be shown to hold when uj = 7 1 = 6 1 / ( 2 + s ) . The solution of eq. fl2^) for 
Q{fl) is unique. The asymptotic value of £ has to be identified with uj^ 1 , and can be used to find 

C(i)=KH. (28) 

(The equation can be solved analytically for one, non physical value of /3 = 3/2.) 

We see then that ro appears as some radius scale in the density distribution of the isothermal 
cylinder. If it can be determined for an actual filament together with a r , we can determine /xo 
through eq.(|2l|) as 

Mo = -^L (29) 
r a o G 

M ° = (2-^W' ( 3 °) 



and, in terms of a± from eq.((2 
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In order to use such relations to determine po for real filaments we have to relate ro to some 
observable property such as the half-mass radius, the projected half mass radius, the radius, i?i/ 2 , 
at which the surface density reaches a certain fraction of its central value, etc.. Take the last, for 
example, if we write the product (2 — (3) 2 rQ appearing in the denominator of eq.(|30[) as qRt/2 we 
have 

4<r 4 



A*0 



qR 1/2 a G 



(31) 



We find numerically that q ~ 1 for (3 = (isotropic case), q ~ 1.5 for (3 = 1/2, and q ~ 2 for (3 = 1 
(radial orbits). 

The second model is contrived to afford analytic solution, but is not less reasonable, in default 
of knowledge of real filaments. It differs from the first in that the velocity dispersions are now 
assumed to increase from the center out as the fourth root of the radius; thus 



(32) 



and f3 is still constant (the projected, central velocity dispersion does not vanish). The constant S 
can be used to construct a quantity 

(33) 



Ga n 



with the dimensions of a line density, which we use to define the dimensionless line density 

C(r) = p{r)/p. 
The Jeans equation can then be written as 



1 djr^p) 13 
p dr r 1/2 



1 djr^p) (3 
p dr r 1/2 



C(r) 



This, in turn, can be written as 



where 5 = (3 + 1/2, and 



din p 
din r 



+ 6 



din p 



din r 



+ 6 



<(r) 



C( r ) = — / p(r')r' dr' . 
A Jo 



Near the origin, where £ ~ 0, we have 



p(r w 0) ar 



(34) 
(35) 

(36) 
(37) 
(38) 



Since, by definition, [3 < 1 the line density always converges at the origin. However, for (3 > 1/2 
the acceleration diverges there, and, in the least, the deep-MOND assumption does not hold. 
These models are thus only good for (3 < 1/2. 



At large distances, where ( w po/p (po, as before, is the total line density) 

p(r — > oo) oc r~( S+T \ 



(39) 
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where r = (//q/ ' (i) 1 ^ 2 ■ As the asymptotic form will be determined by the solution of eqs.(f36|) 
the value of r will be read off the solution, and fiQ will be expressed in terms of /t. These equations 
possess a family of solutions whose members are specified by the choice of the constants, a, 
dictating the behavior near the origin. Defining 



A(r) = a 1 r s p(r) 



and redefining the independent variable 

2-8 



we can write 



x = I — 



87TO 



din A 
din x 



9£(2 - 5) 3 

1/2 

A(x') dx r 



-i 

2-6 



This equation, with the boundary behavior X(x as 0) « 1, has the unique solution 

A(x) = (1 + x 1 / 2 )- 3 . 

Thus, the full density law is 



and 



p(r) 



C(r) = -(2-5) 2 



1 + 



(2-«5)/2 



l -3 



{r/r,) 1 - 5 / 2 
l + ir/r^i 1 - 6 )/ 2 



(40) 

(41) 
(42) 

(43) 

(44) 
(45) 



The radial scale length, r s , characterizing the density law, is arbitrary, and b = ar s 5 is a function 



of r s through eq.(41). In particular we read from eq.pq) that 



M = ^(3-2/3) 2 S 



16 



a„G 



(46) 



As the case is with the first model we can express 5 in terms of different measures of the 
velocity dispersion and scale length of observed filaments. 



For example we can calculate the system average value of a r , or a± 

2 



2 = ,2, = q m Sr^ /2 p{r')dr< 

2-/3 X K r) J ( r r 'p( r ') dr' ' 



The numerator is finite for (3 < 1/2, and we find 

2 



2-/? 



a\ = (a 2 ) = PS^r]' 2 , 



(47) 



(48) 
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P= f-^ U^ ] . (49 ) 



where 

' " (3-2/3) 2 V U1 "3- 2/3, 
The scale radius r s can be easily related, for example, to the (space) half-mass radius, r/,, 
containing half the line density: 

r s = (V2-l) 2 /( 2 -<V (50) 



Thus, in terms of a±, and 



in 2 \ 2 

to = Q-^jt—, (51) 
a Q Gr h 



Q = -^2-l)—?r—^-\ 2 p- 2 . (52) 



where 



4 V ' V2-/3 

For /? = -1/2 P = 3/2, and Q « 6. For /? = P = 40/9^3 ps 2.6, and Q » 2.5. For /? = .25 
P as 4.9, and Q ~ 0.8. For /3 approaching 1/2 Q goes to 0, but then most of the contribution to 
the mean dispersion comes from larger and larger radii. 



4. A preliminary analysis of the Perseus-Pisces segment 



Now we turn to the implications for the Perseus-Pisces super-filament. For reasons explained 
in their paper, ELT actually apply their method to only a segment of the PP bridge of length 
of about 20/i _1 Mpc, and diameter of about 4/i _1 Mpc (h = Hq/WO kms~ 1 Mpc~ 1 ). They 
find a system-integrated, line-of-sight velocity dispersion of a± ~ 430 kms' 1 , estimate the 
B-band luminosity within their chosen section at about 4 x 10 12 /i -2 L Q . Using their a — fi 
relation they get a value of /Uo ~ 8.5 x 10 13 M Q Mpc _1 . All these yield an estimated value of 
(M/L) N ~ 450h(a±/430kms~ 1 ) 2 (M/L) Q . We adopt all the above system-parameter estimates. In 
addition, we need a radius scale, r, for the filament, and choose to normalize rh to 2 Mpc, which 



seems the approximate typical radius from Fig. 1 of Wegner Haynes and Giovanelli (1993) 



We can estimate the MOND M/L value in several ways 
1. Estimate the acceleration a and than correct the Newtonian ELT value by a factor a/a Q . The 



value of a Q as deduced by Begeman, Broeils, and Sanders (1991) , normalized properly to the 
assumed value of Hq is a a ~ 2 x 10~ 8 h 2 cm s~ 2 . The acceleration at r is 

J2 



(7,. 



-1, 



7 



din p 



din r 



and so 



— w 1.5 X 10 



-2, 



430kms 



rh 
2Mpc) 



-l 



(53) 



(54) 



-11- 



Applying this correction factor to the ELT's Newtonian estimate we get the MOND estimate 

E\ ^7 ( a - Y( a ^ Y( rh yV a ° Y l (—\ (5S) 

L J M ~ 7 V430fcms-V \430kms- 1 ) \2Mpc) \2 x 10- 8 h 2 cm s~ 2 ) \LJ @ ' { ' 

2. We can use the above isothermal-cylinder models with all due circumspection. For example, 
from eq.([H]) we can write 

where q is between 1 and 2 for /3 between and 1. With the B-band luminosity estimate of ELT 
this gives 



M\ -.13 -if ^ WvrV yV M ^ 

L/ M ~ 9 \430fcma-V \2Mpc) \2 x 10" 8 /i 2 cm s" 2 / \L/< 
If we use eq. (^) from the second model we obtain an estimator 

/ \ 4 / \ -1 / \ — 1 



(57) 



and 



M0 ~6xl0 n Q -± — 5^ M e Mpc-\ (58) 

pu ^V430fcms-V \2MpcJ \2 x 10- 8 h 2 cm s~ 2 1 

M.\ -W a± Yf^Y'f a ° Y 1 ( MS \ (-<)) 

\430A;ms-V \2Mpc) ^2 x 10- 8 /i 2 cm r 2 j U ' ' °' 



ELT discuss in detail many of the uncertainties that plague the analysis, from doubtful model 
assumptions to difficulties in determining system parameters. We here only stress further some of 
these and add some comments peculiar to the MOND analysis. 

There are several factors that may cause an over-estimation of the velocity dispersion a±: 1. 
Contribution from variations in the mean Hubble velocity along the filament. ELT show how this 
can be largely removed by dividing the filament into segments and taking the dispersion in each 
segment separately; but this was not done for the actual Perseus-Pisces analysis. 2. Clumping 
along the filament that is reckoned without can lead to a substantial over-estimate of a±, as ELT 
demonstrate on model filaments from N-body calculations. The PP segment treated by ELT 
appears knotty to a degree on Fig. 1 of |Wcgner Haynes and Giovanelli (1993)| . 3. motion along 
the filament may contribute if this is not exactly perpendicular to the line of sight everywhere. 
In the case of the PP segment causes 1. and 3. seem to be of little consequence (D. Eisenstein, 
private communication). Such uncertainties in a± are more crucial in MOND than in Newtonian 
mass estimates, because the velocity enters the MOND mass estimate in the fourth power, not the 
second. 



Both the ELT and our method are global, and give the total /i(oo). The luminosity, however, 
is taken by ELT only within a certain projected radius (of about 2 Mpc). The actual filament 
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may well extend farther in radius, and the estimated luminosity, uncertain as it is anyway, may be 
systematically too low because of this, a fact which contributes to an over-estimation of M/L. 

I thank Avi Loeb for bringing the work of ELT to my attention prior to publication, and for 
useful discussions, and Daniel Eisenstein for comments on the manuscript. 
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